Effective-medium theory for infinite-contrast, 2D-periodic, linear composites with strongly 
anisotropic matrix behavior: dilute limit and cross-over behavior 
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The overall behavior of a 2D lattice of voids embedded in an anisotropic matrix is investigated in the limit 
of vanishing porosity /. An effective-medium model (of the Hashin-Shtrikman type) which accounts for elastic 
interactions between neighboring voids, is compared to Fast Fourier Transform numerical solutions and, in the 
limits of infinite anisotropy, to exact results. A cross-over between regular and singular dilute regimes is found, 
driven by a characteristic length which depends on / and on the anisotropy strength. The singular regime, where 
the leading dilute correction to the elastic moduli is an 0{f^^^), is related to strain localization and to change 
in character — from elliptic to hyperbolic — of the governing equations. 



PACS numbers: 46.05.+b,46.15.-x,46.15.Ff 
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I. INTRODUCTION 



Effective-medium approximations (EMAs) for nonlinear 
composite o''^i^i"^'^'^'^i^i^''° (i.e., multi-phase materials), which 
aim to predict their overall (i.e., macroscopic) behavior, are 
pushed to their limits of validity as the nonlinearity and/or 
the heterogeneity contrast become large. ^ Typical examples 
of this sort of phenomenon in continuum mechanics include 
porouS j"''^''^ and rigidly reinforced,'^''* plastic or nonlin- 
early elastic media. In the idealized model of perfect plastic- 
ity, plastic material flow takes place at constant stress intensity 
(the yield stress). In such circumstances, the flow preferen- 
tially concentrates (localizes) in shear bands^^^ Formally, 
these shear bands are closely related to other types of mini- 
mal breakdown manifolds in heterogeneous media (mechan- 
ical systems as well as nonlinear electrical networks) i''^''^i''^ 
However, nonlinear EMAs which address plasticity rely on a 
quasi-equilibrium hypothesis, which means that the charac- 
teristic time of an individual "breakdown" (or slip) event is 
longer than that of wave propagation through the medium^" 
(in nonlinear dielectrics, such conditions are met as well in 
the reversible diode network experiment^'). Plastic deforma- 
tion being a strongly irreversible process, applications of such 
EMAs to plasticity should therefore be restricted to incipient 
deformation in an incremental framework. '^ '^ Nonetheless, 
upon neglecting elasticity and assuming plastic incompress- 
ibility, the quasi-equilibrium hypothesis may extend to full- 
grown deformations, such as in the slip-line theory of perfect 
plasticity, but the governing equations are then hyperbolic.'^ 



Efficient nonlinear EMAs rely on the use of an underlying 
linear comparison medium,^^-^-' which may consist in a "se- 
cant" (isotropic) approximation^^ to the nonlinear response of 
the composite. In the most recent approaches the comparison 
medium is anisotropic, of direction determined by the applied 
field,* and of strength being consistently determined by the 
covariance tensors of the local fields in the phases i^i^i'" How 
these methods cope with localization at the overall level in 



heterogeneous media is not well understood, see Ref. |25| and 
references therein. 

To address this issue, this paper is devoted to the signature 
of incipient localization in an EMA for periodic composites. — 
Because for periodic media efficient methods of solution have 
been developed, ^^-^^'^^'^^-^'''^''^^i^^'^'^ our focus here is on such 
materials. Thus, a system consisting of a two-dimensional 
(2D) periodic lattice of voids embedded in a deformable ma- 
trix is considered. Aimed at understanding the hallmarks 
of localization in the underlying linear medium of nonlinear 
EMAs, we focus on the case of an elastically anisotropic lin- 
ear matrix, of variable anisotropy. ^^•■'^ 

The problem, described in Sec. [Ill admits an exact analyti- 
cal solution in the particular case of infinite anisotropy where 
the governing equations acquire an hyperbolic characteri^^ 
As a consequence, the overall elastic moduli depend on the 
porosity / (the volume fraction of voids) as powers of Z^^^, 
in particular in the dilute limit / — > 0.^^ This result is at odds 
with usual effective-medium results, in which the first cor- 
rection to a homogeneous medium is an 0{f)r^^ due to its 
proportionality to the number of inclusions. 

But for finite anisotropy, the governing equations are ellip- 
tic, and no exact solution is available. The crossover to the 
regime of high anisotropy and its link to localization, of di- 
rect interest for nonlinear EMAs, and more generally for un- 
derstanding the nature of the macroscopic yield transition,'— 
are investigated hereafter. For lack of exact solutions, com- 
parisons are made between: (i) quasi-exact numerical results 
obtained by Fast Fourier Transform (EFT) calculations; (ii) 
outcomes of an EMA for linear periodic media whose predic- 
tive capabilities are assessed; and (iii) the exact results of Ref. 
25. The case of a non-linear (visco-)plastic matrix, of direct 
experimental relevance,^** is examined elsewhere.^'' 

The notation used is as follows: A denotes a tensor of com- 
ponents Aijki ; the sans-serif a is the tensor of components 
(except for the strain and stress e and cr, and the strain po- 
larization T, written in boldface); the boldface a is the vector 
of components a^. A colon denotes a double contraction e.g. 



2 



A : B has components AijmnBmnki, etc. For convenience, 
indices i — x, y or 1, 2 are used indifferently hereafter. 



II. PROBLEM FORMULATION 

The composite, described in Fig. [T] consists of an elastic 
matrix (phase a = 1, of volume fraction c^^^ = 1 — /), con- 
taining a square array of voided cylinders of radius a (phase 
a = 2, of volume fraction c'-^-' = / = ira^). Here and 
henceforth, the size of the unit cell is i = 1. A set of du- 
ality relations'*" allows one to translate the following results 
for the overall behavior of this porous medium, in the context 
of rigidly reinforced composites, which is another interesting 
case of infinite elastic contrast. 

In the composite cr(x) = L(x) : £(x), where L(x) is the 
position-dependent elasticity tensor, of components Lijki — 
Lijik — Lkiij- The equilibrium equation c)iCTy = holds, and 
the strain derives from the displacement u as = {diUj + 
djUi)/2 (small perturbations are assumed). In two dimen- 
sions, ea;a; = d^Ux:, £yy = OyUy, SLUd E^y = {O^Uy + SyU^) / 2 . 

In the voids, L(x) = L^^^ = 0, the stress vanishes, and the 
strain is arbitrary: any continuation matching the displace- 
ments at the voids boundaries is admissible. Only its volume 
average over the void is physically relevant. 

The matrix material can be thought of as a "mixture" of two 
basic types of anisotropic media;-^ (i) one where the eigendi- 
rections of anisotropy coincide with the reference axes of unit 
vectors ei = and 62 = e^; and (ii) one where they coin- 
cide with the diagonals (see Fig.[TJ. Such a medium is invari- 
ant under the dihedral point-symmetry group D/^^ Then, its 
elastic tensor L^^^ is of the form 
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2Li2l2^^^ + {Lllll-Lii22W\ (1) 



PS 



where J, E^^'^^ are mutually orthogonal projectors defined by 
(I, of components 5ij, is the 2x2 identity matrix): 



J = (1/2) I® I, 
SS'PS = (1/2) 



SS,PS^gSS.PS^ 



(2a) 
(2b) 



The identity is I = J + E^^ + E^^. The definitions involve the 
so-called simple shear (SS) and pure shear (PS) eigenmodes 
of deformation: 



cc 

e = ei (g) e2 + 62 (8) ei, 

pc 

e = ei (g) Gi - 82 (g) e2, 



(3a) 
(3b) 



such that E^^'^^ : e^^'^^ = e^^'^^. Their eigenvectors are re- 
lated by a 45° rotation (see Fig. [T]i. Also, J : I = I. This 
decomposition relates to that of a 2 x 2 symmetric tensor a 
into one equibiaxial (i.e. compressive) mode and two orthog- 
onal shear modes: 



a = a„ I + ass + aps 



(4) 



of respective components = {clxx + o-yy)l% «ss = O'xy, 
aps = (a-xx — aj/i/)/2- Thus, in the matrix we write: 



.(x) = L(i) = 2k J + 2A E'^'^ + 2^ E- 



K is the bulk compressibility modulus, and A, [i. are in-plane 
anisotropic shear moduli. With this medium of a special or- 
thotropic type, the interpretation of the 2D problem as a lim- 
iting one of plane stress (cr^rz = <^yz = ^zz — 0, £zz 7^ 0) or 
of plane strain (Exz = ^yz ~ ^zz = 0, (Tzz 7^ 0) is irrelevant 
from a theoretical standpoint, though the expression in terms 
of K and of fj, of the Young modulus and Poisson ratio relative 
to the pure shear mode differ in both cases. 

With applications to volume-preserving plastic deformation 
in mind, this study mostly focuses on the limiting case of an 
incompressible matrix for which k — 00. Introducing in this 
limit the dimensionless anisotropy ratio k = A//i, the medium 
is isotropic when fc = 1, and is infinitely anisotropic when 
either fc = or fc 00. In each of the latter limits, the 
medium possesses one infinitely hard, and one infinitely soft 
eigenmodes of strain: when fc = (i.e. X = or fj, — 00) 
the medium is soft for SS loadings, and resists PS loadings, 
whereas when k = 00 (i.e. A = 00 or /i = 0) the medium 
is soft for PS loadings, and resists SS loadings. We accord- 
ingly call these loading modes "hard" and "soft" hereafter. 
This model provides a convenient framework for studying the 
coupling between porosity and localization. 
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FIG. 1: Left, periodic porous medium with unit cell and reference 
axes. Right, black arrows depict eigenmodes of strain: simple shear 
(SS), and pure shear (PS). 



Volume averages over the whole system are denoted by 
brackets (•). Hereafter, (•)'■"■' stands for a volume average 
over phase a. By definition, the effective (overall) tensor of 
elastic moduli, L, relates the macroscopic strain e = (e) to 
the macroscropic stress a- = (cr): 



L : e. 



(6) 



PS 



(5) 



The square void lattice also being invariant under D4, L is of 
a form analogous to (|5]l, where k. A, /i, are replaced by the 
effective moduli k, A, Jl. The latter are the main quantities of 
interest. They depend on k. A, /i, and /. Even when k = 00, 
the effective modulus k is finite for the porous medium. Then, 
the normalized moduli X/X,Jl/fj, depend only on k, and on /. 
Convenient normalizations for k are k/ii when A 00, or 
k/X when /i — > 00. 



III. FULL-FIELD FFT APPROACH 

A. Numerical method 

Full-field numerical solutions of the problem are ob- 
tained using the Fourier transform method,^"* applied to linear 
composites. The method amounts to solving iteratively the 
Lippmann-Schwinger equation for the strain^^ 



£(x) = £ + y" d22/G(x-y) :T(y), (7 
t(x) = [l(x) -L^")] : £(x), 



a) 
(7b) 



where L^^^ is some arbitrary background elastic tensor. The 
position-dependent elastic tensor of the medium, L(x), is 
( — L^^^^) in the void and L*^^^ in the matrix. In all the nu- 
merical calculations of the paper, the latter is assumed nearly 
incompressible with k ~ lO'^, and no appreciable differ- 
ences were observed for n ~ 10^. The tensor G is the pe- 
riodic Green function of the background medium, such that 
J d^xG(x) = 0. In the method, the convolution in Eq. ( iTat 
is evaluated in Fourier space, whereas (iTbl i is computed in di- 
rect space. The system is finely discretized as a L x L array 
of pixels. The bad iterative properties of O are alleviated 
through various improvements allowing for high or even in- 
finite contrast.^^'^*' "*^ These schemes are used here. Fast con- 
vergence is achieved by taking L'°) of the type ([5]l, with the 
same anisotropy ratio k as L^^\ but with considerably lower 
moduH, namely ~ 5. 10"'* and ~ A"/A ~ 0.2 

(not necessarily optimal values). The Fourier transform of G 
reads^ 



y •'J sym 

where { jsym indicates a symmetrization so that Gijki 



(8) 



G 



klij 



— Gjiki, and where (q) = qk qi is the acoustic ten- 
sor. 

Calculations are carried out for various anisotropy ratios 
10""* < k < 10'', and porosities < f < fc using FFT 
routines. Three sizes L = 512, 1024, 2048 are considered to 
monitor size effects. The smallest one leads to results with 
satisfactory convergence properties, except in cases of high 
anisotropy where a better resolution was required to render the 
fine structure of the field patterns. We used L = 2048 when- 
ever an appreciable difference was found between i = 512 
and 1024. 

Once the strain Sij (q) is computed, the displacement u(q) 
is deduced from (q ^ 0)4^ 

Uxiq) = -i {Qx [exxiq) - Syyiq)] + 2qye^y{(i)} /g^ 
Uy{q) = i{qy[e^:,r,{q) - eyy{q)]-2q^e^y{q)} /q'^ 

Only SS or PS macroscopic strain loadings are considered 
(sm = 0). Other shear states follow from linearity. For 
both modes, the linear elastic problem is solved for various 
anisotropy ratios < fc < oo. Effective moduli are computed 
using one component at a time, e.g. A = (cTxy) /[2{exy)]- 
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FIG. 2: Structure of unit-cell field patterns for high anisotropy. Left, 
pattern for SS loading and k — oo. Right, pattern for PS loading 
and fc = 0. In these figures, v = void; b, c, d — deformation bands 
in the matrix (not intersecting, intersecting far from the void, and 
intersecting around the void, respectively); a — remaining parts of 
the matrix. 



B. Overview: displacement and stress maps 

Typical displacement and stress maps obtained by (iso- 
choric) FFT calculations are as follows. Since first-order in- 
finitesimal displacements are used, it should be borne in mind 
that however singular, the displacement patterns are at most 
incipient ones. In limits of infinite anisotropy fc ^ 0, oo, the 
structure of the solutions tends towards that schematized in 
Fig. |2l being organized into bands of width one void diame- 
ter. Two remarkable types of structures are found, depending 
on the loading direction and on the type of anisotropy. They 
differ essentially by the presence of zones in the matrix where 
the bands cross (denoted by c in Fig.|2]i. 

With the above mentioned caveat, the following approxi- 
mate symmetry holds between the maps:^^ 



7?.45o(void lattice) <^ 



fc ^ 

PS loading 



1/fc 

> SS loading, 



(9) 



where the TZ symbol denotes a 45° rotation of the lattice of 
voids, with all other parameters (material constitutive law and 
loading) conserved. These field structures, already revealed 
by the analytical calculations of Ref.|25|for infinite anisotropy 
(to which we refer the reader for further details), are retrieved 
here for finite, but high, anisotropy. 

Table U displays full-field calculations of the reduced (peri- 
odic) displacement field u*(x) = u(x) — e • x, indicated by 
arrows, superimposed on a representation of the unit cell de- 
formed using a rescaled displacement f3 u* (x), for anisotropy 
ratios fc ~ 10^^, 1, 10'^ in SS and PS loadings, for a moderate 
porosity / — 0.1. To highlight the deformation pattern, the 
magnification factor f3 lies between 1 and 10. Lighter grey 
tones in the deformed matrix indicate regions subjected to a 
strong extension. The unit cell is repUcated in order to empha- 
size the displacement ("flow") pattern. For low and high fc, the 
features of the displacement maps are in agreement with the 
exact results derived at infinite anisotropy in Ref. |25l which 
they enlighten. 

The flow pattern is organized in closed convection cells of 
square shape, delimited by black boxes. Two types of cells, 
rotated 45° with respect to the Cartesian axes, and related 
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TABLE I: Reduced displacement field u* (arrows), and resulting elastic deformation of the unit cell (to lowest order of perturbations), in SS 
and PS loadings for anisotropy ratios k — 10"'', 1, 10''. An enlargement of a void is shown in map (D). 



k = IQ-^ k = l fc = 10^ 




by a mirror symmetry, suffice to account for the flow pat- 
tern in SS (maps A, B, C). As a consequence, and due to the 
high anisotropy, the edges of the unit cells in (A) and (C) un- 
dergo non-zero and quasi piecewise-linear deformation. On 
the other hand, four types of convection cells, aligned along 
the Cartesian axes, related by mirror symmetries with respect 
to these axes, and fully enclosed within one unit cell, are re- 
quired to produce flie flow pattern in PS (D, E, F). 

Compared to the k = 1 isotropic solutions of (B) and (E), 
solutions for highly anisotropic situations are either: (i) lo- 
calized in strain, with a displacement field discontinuous at 
places [maps (C) and (D)]; (ii) localized in stress, with contin- 
uous displacement as in (A) and (F). Strain-localization arises 
whenever loading along a "hard" mode takes place. Then, 
the highly anisotropic medium resists most the applied strain 
and undergoes both a high induced stress and a weak induced 
strain. In the limit of infinite anisotropy, a rigid "block slid- 
ing" incipient pattern results, where the flow is organized in 
bands of width one void diameter (see also Fig. |2]i, where 
the tangential component of u* is discontinuous, and where 
strain concentrates as Dirac distributions along the sliding 
lines. This pattern is tantamount to a breakdown mechanism. 
In turn, block sliding leaves four incipient voids in the matrix 
in (C) and (D), at locations where the sliding lines intersect at 
90°. One such void is enlarged in (D). One important differ- 
ence between cases (C) and (D) is that in (C), the flow bands 



(of width one void diameter) cross inside the matrix due to 
their 45° orientation. Flow redistribution then takes place in 
the intersection zones. On the contrary, in (D) such zones do 
not exist in the matrix, and flow redistribution requires a non- 
zero displacement component normal to the band boundaries. 
As a result, the gradient of the tangential component of u* is 
higher in (C) than in (D). 

Solutions with a continuous displacement field are ob- 
tained instead when loading is applied along the soft defor- 
mation mode. The remaining deformation mode being harder, 
this leaves less possibilities for easy deformation than in the 
isotropic case; this explains why the unit cells of (A) and (F) 
are much less deformed - a magnification /3 = 10 is used - 
than the cells (B) and (E) of the isotropic material - plotted 
with f3 = I. 

TableHJ displays, for / = 0.1 and for increasing anisotropy 
ratios k, maps of the independent stress components am, crp 
and (j±, under SS and PS loadings. Hereafter, the maps are re- 
ferred to by their individual number (1 to 30). Each map goes 
along with its own field scale at its right, in correspondence 
with the color scale at the extreme right of the rows. The 
"parallel" (||) and "perpendicular" (_L) notations refer to the 
"direction" of the applied macroscopic loading. The SS and 
PS shear components of the stress are defined in Eq. (01). In PS 
loading, we have (T[| = a^^, a± = a^^, whereas in SS load- 
ing: (III = (7^^, (j± = a^^ (herafter, a similar notation is used 
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TABLE II: (Color online) Parallel (o-y), transverse (a±), and mean (am) stress field maps for SS and PS loadings, with anisotropy ratios 
k = 0.01, 0.2, 1, 5, 100 (porosity / — 0.1). The stress fields are rescaled such that (cr||) — 1 

fe = 0.01 fe = 0.2 k=l k = 5 fe = 100 



Simple shear 




-3.14 -2.73 -2.87 -3.10 - ■ « «« 



for strain components). In both cases, the volume average of 
the non-parallel components vanish: {a±) — {(Jm) = 0. The 
maps display rescaled stresses, such that ((T||) = 1. Due to 
linearity, the strain fields are the same, up to a change of scale 
(although the scales are different in the parallel and perpen- 
dicular directions due to the anisotropy). 

The following observations are relevant to the regime of 
high anisotropy, where the stress patterns follow that of Fig.|2] 
The zones where bands cross depicted in this figure [either 
in the matrix (zones c) or close to the voids (zones d+v)]. 



are places of additive screening or enhancement of the stress. 
Thus, the parallel stress in zone c of map 20 reaches its highest 
values there, and is twice that in the two crossing bands (how- 
ever, a much higher transverse stress is encountered in the im- 
mediate vicinity of the void, see map 25). In a similar way, the 
vanishing stress in zone c in map 5 is the difference between 
the stresses in the bands. Two remarks, strictly valid for infi- 
nite anisotropy, are in order at this point: first, zones of van- 
ishing stress are squares, of size determined by the void cross 
section transverse to the bands, so that the disk-like shape of 
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the voids is no longer relevant; second, the build-up of zones 
of zero stresses (i.e., analogous to porous zones) in the matrix 
in SS loading leads to an effective doubling of the porosity in 
the effective shear modulus A at infinite anisotropy, whereby 
an effective "close packing" threshold, twice as small as the 
geometric one,^ is reached as / increases, leading to a "me- 
chanically advanced" percolative behavior. As a consequence, 
A decays rapidly with /, see next section. 

The stress is less singular than the strain in the limit of in- 
finite anisotropy. Indeed, in a strain-localized situation (load- 
ing along a "hard" mode [maps 5, 16]), the displacement is 
discontinuous. Accordingly, the transverse strain has Dirac 
singular components along the band frontiers. They abruptly 
change sign at the special points (±a, 0) and (0, ±a) on the 
void boundary in PS and at points (±a, ±a) / \/2 in SS, where 
a is the void radius. Because of the stress-strain proportional- 
ity, these strain singularities can be traced in maps 10 and 21. 
However, since the perpendicular stress vanishes in the limits 
k ^ 0, oo, so do its Dirac singularities, as shown by the small 
values on the scales. The special points, termed hot spots in 
Ref. 25 are points of extreme matter separation, or crushing, 
which bear the main cost of the "block sliding" patterns. On 
the other hand, the incipient secondary voids in maps (C) and 
(D) of Tab. U appear (somehow paradoxically) as regions of 
moderate stress levels. 

More generally, the stress field undergoes the following 
types of singular behavior in the limiting cases of infinite 
anisotropy: 

— loading along a hard mode: discontinuous a\\ component 
along band frontiers in the direction normal to the frontiers, 
with finite jump, accompanied by hot spots at the void bound- 
ary [maps 5, 16]; 

— loading along a soft mode: discontinuous derivative of a\\ 
in the same direction, with infinite jump [maps 1, 20], and 
discontinuous a±^ with infinite jump across the band frontiers 
(case of loading along a soft mode) [maps 6, 25]. 

The mean stress is always singular with the most singular be- 
havior: it has the singularity of the parallel stress in the case of 
loading along a soft mode, and the singularity of the transverse 
stress in the case of loading along a hard mode (but the mean 
strain vanishes in the limit of an incompressible medium). 

IV. ANALYTICAL EFFECTIVE MEDIUM APPROACH 

Nemat-Nasser proposed'^ an approximate (dipolar) Fou- 
rier-mode approach to the periodic problem, which proved 
excellent for isotropic componentSi^^iS We apply it to the 
anisotropic case. More accurate schemes going beyond the 
dipolar level, however less suitable to analytical treatment, are 
available.2i22i32 

The approach is as follows. Consider first the general case 
of a binary composite of volume V — > oo, the inclusions of 
which have an elastic tensor L*^^', and set (5L = L^^^ — L*^^^ 
The characteristic function Xoo of an infinite periodic array of 
identical inclusions, of characteristic function x, is Xoo (x) = 
x(x ^ Ti), where are lattice vectors. Then, L = L*^^' + 
Xooi^L- Equations Q, (O apply, with L(°) = L^^), and r = 



XooSh : e. Multiplying ( fTal l by x, integrating over V, and 
assuming homogeneous deformation in the inclusions so that 
XooE ~ Xoo (e) results in an expression of (e)^^^ in terms 
of the Hill depolarization tensor^^'^^ of the lattice 

^ = ^y^) J d'a:dVxoo(x)G(x-x')Xoo(x') 
= -^^/(§2E^"-^G(q)|(e--)/ (10) 
= -c(2) 5]'G(q)|(e-'i)/, (11) 

q— 27rp 
pGR.L. 

where (•)/ denotes a volume average over one individual in- 
clusion of volume Vi. The sum in ( fTol ) is over lattice sites. The 
last equality stems from the Poisson summation formulai^ 
The primed sum is over nonzero reciprocal lattice vectors 
p ~ {px,Py) (with integer components). Carrying over the 
obtained (e) to the volume average (L : e) computed from 
(l7]l, and using definition ([6]), entails the effective elastic tensor 

L = L(^) +c(2)(5L : (I + P : (5L)"\ (12) 

The formula for the void lattice with L'^^ — follows. 

Equ. (fT2l i is of the Hashin-Shtrikman (HS) variational 
|.ypg^55j56j57 particular, the void lattice version provides an 
upper bound to the exact result.^" Interpreted in the frame- 
work of multiple-scattering theory,^^'^^ Equ. ( [T2l i states that, 
at the dipolar level, the effect of the inclusion lattice on 
the homogeneous matrix amounts to that of non-interacting 
identical "equivalent" inclusions in proportion c^^\ charac- 
terized by a T-matrix provided by the term following c*^^^ in 
(jl2j|.,58j59 furthermore, each of these equivalent inclusions pos- 
sesses a lattice, c*^^^ -dependent, microstructure, the influence 
of which is lumped in P. To make contact with other types of 
effective-medium expressions, separate P into one and two- 
body contributions^** by writing P = Pi — c'^'P2, where Pi is 
the r = term in ( fTOl i. A similar decomposition is proposed 
in Ref.i60|in the context or random composites. Then, intro- 
ducing = L — L^^^ (fT2l i takes on the familiar HS form 
with pair correlations"*^ T = c'^^^Ti where, 

T EE SZ : (1 + P2 ■■ sVj \ Ti = 5L : (I + Pi : ^L)"\ 

(13) 

This decomposition proves useful whenever some eigenvalue 
of Ti should blow up. Then, the corresponding eigenvalue of 
Sh is simply provided by that of —^2^^ above expres- 

sion makes clear. 

For cylindrical voids of radius a, with Ji the Bessel func- 
tion, (e^'i'')/ =2Ji{aq) /{aq). Setting 

m = £=X/k, (14) 

one finds from ([Hi, (El that Pi^ki = {QipqiRjpqkjsym, where 
M = JJ - (1 + m)E^^ - (1 + ^)E^^ where "sym" denotes 
a symmetrization with respect to indices {i, j) and (fc, I), and 
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where 



TT 5Z 



pe 



(27rap) 



(15a) 



A(p) = A(l+m)(p2_p2) +4Ai(l + £)p2p2. (15b) 

The reciprocal lattice is a square lattice. Hence Q is also in- 
variant under D4. Being completely symmetric, it is of type 
([U with L1122 = ii2i2 and is determined by two independent 
scalar lattice sums only. One obtains: 



1 



+ i) 



where, after having reduced the lattice sums to sums over the 
positive quadrant. 



5a 



Px>0 P 

Py>l 



{pI-pVi 



ApIpI 



These sums bring in the anisotropy parameter: 

A:^[(l + m)A]/[(l+^)/i] 



(17) 



(18) 



which reduces to A//i in the incompressible limit k ^ 00. We 
remark in passing that 



kS\ + S"^ = S2{a) 



J2 [Ji{'27:ap)/pf (19) 

Px>0,Pa>l 



is independent of A;. After some algebra, one arrives at 
1 



^ , (l + m)5A]ESS 
(20) 



The one-body Pi is read from this expression, provided that 
S\,^ are computed in the continuum limit, by making the sub- 
stitutions J2 ~^ iliine_>o d\/{2Tr)'^, p ^ q/(27r) in 
([TTT i. Then (in the continuum limit), S'a,^ 



S 



111 



l + \/fc' 



(21) 



Eqs. ( fTTT i show that 5a blows up when k ^ due to the 
contribution of the Cartesian axis px — 0. On the other hand, 
Sfj. remains finite or goes to zero in all cases. 

The limit of an isotropic matrix where k = 1, £ ~ m ~ 
K/fi provides Pi = [2r7iJJ + (1 + 2m)K]/[4^(l + to)], where 
K = E^^ + E''^ . This expression can be recovered directly 
from ([8]) and from the usual definition in terms of an angular 
integral^' Pi = - J d^qG{q) / {27^) where q = q/q (the 
independence wrt. x stems from the rotational symmetry of 
the voids). 

From ST% . ( I20I ), the effective moduli of the void lattice 



[{Sx + Sf,)I + S^E^s + SxE^^] , (16) read, with / 



42). 



k/k = l-f/{l-[{X/^l)Sx + S^]/ {!+£)}, (22a) 
A/A = l-//{l-A:[S'A + m5^/(l + TO)]}, (22b) 
m/m = 1 - // {1 - [S^ + mSx/{l + £)]} . (22c) 

Henceforth, incompressibility is assumed for simplicity so 
that fc = A//i from now on, unless explicitly stated. 



V. RESULTS 
A. Effective moduli 

The numerical results at various values of k and / discussed 
in this section are obtained using brute force numerical com- 
putations of the sums Sx.fj., with convergence checks. The 
sums Sx,fj, are conditionally (and slowly) convergent and the 
following suitable prescription is used. Sums are carried out 
over concentric square shells of points 5„ ~ {{px,n) \ < 
Px < n-l}U{{n,py)\ 1 < py < n} for 1 < n < N, with iV 
is sufficiently large. Huge numbers of terms are required for 
accuracy, especially in the dilute limit. 

Figs. |3] show comparisons between the effective moduli 
computed numerically from the above maps, and the EMA of 
Sec. lIVK indicated as HS in the plots), for / = 0.1. The agree- 
ment is excellent near the case of an isotropic matrix k — 1 
(as is expected for such a small porosity), but also up to high 
anisotropy. In all cases, the EMA is seen to provide an upper 
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ratio. The point k^,, as estimated by the EMA, is represented 
vs. / in Fig.|4](no attempt has been made to use full-field cal- 
culations for computational cost reasons). The curve 
shows that crossing only occurs for porosities / < ~ 0.13: 
as / increases, the curve A/ A in Fig.[3la) goes down to zero 
faster than Jl/ fi, while the crossing point shifts to the left un- 
til it vanishes. Remark that fc* < 1 whenever it exists. For 
A: < fc* < 1, the matrix is more anisotropic than the compos- 
ite; the inverse situation prevails for fc > fc*, and in particular 
for / > /* where fc* = 0, so that void-induced anisotropy 
dominates in this regime. 



FIG. 4: Anisotropy ratios fci and k, vs. porosity / in the HS-EMA. 
Incompressible matrix. 



bound for the corresponding full-field estimates (a property of 
the Hashin-Shtrikman approach). 

Near fc = 1, the orientation of the void lattice makes the 
medium harder under PS loading, than under SS loading [i.e., 
//(/c = 1) > A(fc — 1)]. Indeed, the anisotropic matrix can be 
thought of as containing rigidifying fibers (of strength /i), ori- 
ented at 45° along the diagonals, that resist PS deformation, 
and fibers (of strength A) oriented along the Cartesian axes, 
that resist SS deformation. In the void lattice, the nearest 
neighboring voids, and consequently the largest "directional 
damage", lie along the Cartesian axes, which explains the dif- 
ference. We emphasize that while this observation remains 
true as fc ^ 0, the situation changes as fc — > oo: in this 
limit, due to the presence of the fictitious voids produced by 
band crossing, mentioned in the previous section, the nearest- 
neighboring "voids" become located along the diagonals, so 
that the PS direction becomes, for k higher than some value 
k = ki (discussed below), the most damaged one, hence the 
softest. 

As fc — > (Fig. [3^), the curvature of the plots indicates that 
Jl increases slower than /i, whereas A decreases slower than 
A. As fc oo (Fig. [3]3), the parts played by A and fi are 
reversed. Moreover, the compressibility modulus k decreases 
as anisotropy increases, in a way comparable to the hardest 
shear modulus (Figs. [3]; and d). 

At high anisotropy fc, 1/fc < 0.1, discrepancies between 
full-field calculations and the EMAjirise for the hardest shear 
modulus (i.e., Ji when fc 0, and A when fc oo), whereas 
the softest one remains extremely well reproduced. This may 
indicate that the lattice sums have problems dealing accurately 
with the effect of second-nearest-neighboring voids. Indeed, 
the softest direction is always the one where the voids (real, 
or fictitious) are nearest-neighbors, whereas the hardest one 
corresponds to second-nearest-neighbors. Note that k, which 
behaves as the hardest effective modulus, suffers similar dis- 
crepancies at high anisotropy. 

To discuss the crossing of the curves that takes place in 
Fig. [3^, consider the effective anisotropy ratio fc = X/Jl. 
Crossing occurs when fc(fc*) = fc*, for some fc = fc* (/) where 
the overall medium and the matrix have the same anisotropy 



The other remarkable anisotropy ratio is the aforemen- 
tioned fci, defined by the equation fc(fci) = 1, where the over- 
all behavior is isotropic in the plane. This point, also repre- 
sented on Fig. m exists at least up to high porosity values. 
However, since the EMA is expected to fail around / = 0.5 
(see below), the irrelevant part of the fci (/) curve is sketched 
with dashed lines in Fig.|4] In the relevant porosity range, the 
fact that fci (/) > 1 indicates that the matrix needs to be made 
harder along the SS (A) mode than along the PS (/i) mode in 
order to reach isotropy, so as to compensate for higher soften- 
ing in this direction due to newly appearing nearest neighbor- 
ing voids, as is explained above. 



Fig Sillustrates the behavior of the moduli with the poros- 
ity /, for finite anisotropy ratios fc = 0.01, 1, 100, together 
with the exact results of Ref. 'zs' at fc = 0, oo. The exact 
curve for at fc = oo, almost superimposed with the EMA 
curve for fc = 100 in (b), is available up to / — tt/S only,— 
Firstly, the EMA is again seen to systematically overesti- 
mate the moduli. Next, all the elastic moduli must vanish at 
least at the geometrical close-packing threshold of the voids, — 
f ^ fc — tt/A ~ 0.78, and possibly before. The EFT points 
in Figs.|5^ andlS]? are consistent with this fact, whereas the 
EMA fails by producing non-zero results at this point. This is 
not surprising, since EMAs of the HS type are known not to 
be able to account for percolative-type behaviori^ Moreover, 
the exact result in Fig.|5^ for fc = oo shows the shear modu- 
lus in the hard direction. A, to vanish at / = /c/2 due to the 
fictitious voids produced by band crossing. Accordingly, for 
large but finite fc, A decreases rapidly with / up to / = /c/2, 
then with a lower slope up to / = /c. The EMA again fails to 
account for the threshold at /c/2, although the local minimum 
of A at / ~ 0.55 in Fig. |5^ may indicate that at least part of 
the phenomenon is captured by the dipolar lattice sums. Inter- 
estingly enough, when available, the exact results for infinite 
anisotropy at fc = (resp. fc = oo) are seen to provide tight 
lower (resp. upper) bounds to the effective moduli for all val- 
ues of fc, and in particular to the isotropic case fc = 1. As 
far as the effective moduli are concerned. Fig. |5] clearly shows 
that the EMA can be trusted quantitatively up to / = 0.30 at 
most, and is qualitatively reasonable (as long as the matrix is 
not too anisotropic) up to / = 0.5. 
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FIG. 5: Effective shear (A, /i) and compressibility (k) moduli vs. porosity / for various anisotropy ratio k. Incompressible matrix. Comparisons 
between the EMA for = 1 (solid), and k = 0.01, 100 (dotted); FFT results for k = 0.01 (triangle dots), k = 0.01 (circle dots), and k = 100 
(losange dots); and exact results at fc = 0, oo taken from Ref.l25l (dash-dotted). 



B. Continuous transition in tlie dilute limit f <^ 1 

1. Finite anisotropy 

For a finite anisotropy ratio k, the dilute expressions for the 
effective shear moduH at sufficiently small / are read from 
expressions ( I22bb . (I22cb with Sa,^ replaced by the one-body 
contributions 5ia and Si^ defined in ( |2TI ). For the incom- 
pressible medium, the HS estimates of the shear moduli are: 

A/A = l-/(l + VI) + 0(/2), (23a) 
m/m = 1 - /(I + + O {f) . (23b) 

As to the effective compressibility modulus, the incompress- 
ible limit leads to the situation described below Equ. ( fT3] ), 
where one eigenvalue of Ti blows up. This requires us to go 
beyond the one-body approximation. However, Equ. (lA.Sbb in 
the Appendix n shows that 52(a) in ( fT9l ) is exactly S2 — 1 — f 
for / < 7r/4. Replacing, e.g. 5^ by 6*1^ + 0{f) and com- 
puting Sx via (fT9] l, then letting k ^ 00 in ( I22ab [with k read 
from ( fTsT ll provides: 

K= \/V// + 0(l). (23c) 

Remark that "extended" dilute approximations, which extrap- 
olate the above formulas for moderate anisotropy to finite (but 
small) porosities, result from taking P2 = Pi in ( fT3l ), i.e. from 
using in ( |22] |: 

5^,a = 5^"a = (1-./)5i^,a. (24) 

This amounts to assuming pair correlations between the voids 
dictated by the void shape,*" and provides HS formulas of the 
"classical" type in which the lattice structure is ignored. 




(«] (bl (c) 



FIG. 6: (Color online) FFT computations of eps — e\\ in PS load- 
ing, for a material with strong anisotropy ratio k = 10""^ ^ 1. 
Porosities from left to right: / — fr/W (regular "dilute" regime); 
/ = fr (crossover regime); / = 10/r ("strongly anisotropic" 
regime), where /,. ~ fc/7r is the cross-over porosity. Incompress- 
ible matrix. Black represents the highest field values (arbitrary color 
scale). 



2. Infinite anisotropy 

After the incompressible limit k ^ 00 is taken, the limit 
of infinite anisotropy A; ^ is obtained by letting fi ^ 00 in 
A, and A ^ in /i. Conversely, k 00 requires ^ in 
A, and A ^ 00 in /i. In these limits, the sums ^a, 5^ in ( fTTl ) 
are computed in the Appendix„ The obtained shear moduli A, 
Jl are compared to the exact results of Ref. One finds for 
k^O: 

(25a) 
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while exact expressions are: 



3/2 



6 8 



^ = l-i'l 
For k 



1/2 



(25d) 
(25e) 



oo, the HS estimates read: 
3^2 



^ = 1- 



16\/2 V^r 



/ 

512 Vtt 



K 

A 



16%/2 /vr 
[j 



1/2 



- 2. 



fO(/3/2), (26a) 
256\ //^ ' 

+ 0(/^/2), (26b) 
(26c) 



whereas exact expressions are: 



A 
A 

M 



1-/ 



1/2 



32 // 



3/2 



1 



(26d) 



4 ) /'■(26e) 



The above comparisons show that the HS estimates do an ex- 
cellent job of capturing the presence of half-integers powers 
of / in limits of infinite anisotropy at lowest orders in the di- 
lute limit. Moreover, even when the numerical coefficients are 
not exact, they are close to the exact values. The less singular 
character of A in (I25ab . (I25dl l when fc — + [resp. Jl in (I26bb . 
(I26eb when k 00] is discussed in Ref. 251- 



can be written with a summand proportional to [l+i4'7(p)] ^, 
where the dimensionless quantity 7(p) reads 



7(P) 



^pIpI 



(27) 



Singling out the contribution of the main diagonal to S^, 
the remainder of this sum can be brought down to a sum 
over Pa; > 2 and 1 < Py < Px — 1, in which 1/4 < 
7(p) < (Pk/I)^- Hence, 7(p) ^ p^ so that K^^p) pro- 
vides an appreciable fc-dependent contribution only for p > 
I/VK ^ 1/^/k. Moreover, [2Ji{x)/x]'^ is appreciable only 
when X < 2. In terms of p, this reads p < l/(7ra), see 
(fTTb . Hence fc-dependent terms contribute only provided that 
1/Vk < p < l/(7ra). In turn, this is possible only if 
^{^) S 1- For ^(fc) > 1, a fc-independent regime instead 
takes place in 5^. 

From a physical standpoint, the length scale ^ represents 
an effective inclusion size. Fig.|6]indeed displays three maps 
of the parallel strain field in PS loading, computed by FFT 
at fixed anisotropy ratio fc ~ lO^'^ with varying porosity 
/ ~ /r/10, fr and lO/r- It is seen that localized shear 
bands develop from the void as porosity increases. At regime 
change, they coalesce and span the entire medium. The void 
can be considered as an isolated inclusion only for f < fr- A 
similar effect takes place for high fc values. We checked nu- 
merically that in both cases, before coalescence, the strain in- 
tensity in the bands decays exponentially as e cx exp ( — 6 r/^) , 
where r is the distance from the void, and where & is a numer- 
ical coefficient of order one. 



C. Average fields and standard deviations 

1. General considerations 



3. The dilute transition 

Obviously, a cross-over takes place between sets (l23T l on 
the one hand, and (|25l l, (|26] l on the other hand. Balancing 
the "extended dilute" sum S'^f ^ with S^"^^" [equ. ( TObj il, 
then with [equ. ( |A.9dH and solving for fc, provides a 

discontinuous cross-over porosity fr{k) curve which defines 
in the (/, fc) plane boundary lines between dilute and high- 
anisotropy regions. Owing to the approximations at play, 
this boundary cannot be trusted for fc of order one (for this 
reason we do not display the curves). On the other hand, 
we find fr{k) ~ (97r3/lG24)fc ~ k/n for fc < 1 and 
frik) ~ (97r3/512)fc-i ~ 2/(fc7r) for fc > 1. 

Due to the relation / — tto^, the cross-over porosity in the 
highly anisotropic regime stems from a length scale ^(fc) such 
that S. ^ ajk^l"^ for fc < 1, and i ~ afc^/^ for fc > 1. From a 
mathematical standpoint, these length scales originate from a 
scaling property of the lattice sums. We focus here on the case 
fc ^ 0. The case fc ^ 00 can be discussed by adapting this 
argument. Introducing K = fc/(l — fc), the sum 5*^ in ([TtT i 



The first two moments of the fields are required for appli- 
cations to non-linear EMAs, and can be consistently com- 
puted from any linear homogenization estimate.^ Hereafter, 
ei"'' = (e)*'"'' /e denotes the phase average of a strain com- 
ponent e, normalized by the applied macroscopic field. Like- 
wise, we denote by SD^^^e) its standard deviation (SD) in 
phase a , normalized by e. Similar notations apply to stress 
components. 

The phase-averaged fields in the porous composite are de- 
duced from the set of equations 

l-.e = (1 -/)L(i) : (£)(i\ (28a) 
e - (1 -/)(£)(!)+ /(e)(2). (28b) 

Moreover, assuming single mode-loading, the second mo- 
ments in each phase are obtained by taking a derivative of the 
strain energy with respect to the elastic moduh of the phases, 
asi 

(4,ss,ps)^"^ = ^^|^(em,ss,ps)^ (29) 
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where Lis k (resp. A, /i) when the index in the l.h.s. is m (resp. 
SS, PS) and where is (resp. A^"), ^(")) when the 
index the r.h.s. is ni (resp. SS, PS). The variances follow. If 
need be, the incompressibility limit is taken after these quan- 
tities are computed. 

TableHnidisplavs for f = 0.1 the normalized phase-average 

—(a) 

Strains Se for a = 1, 2 and SDs of the strain and stress 
components in the matrix, as computed by the EMA and by 
full-field calculations. The overall agreement is again excel- 
lent, the most important observed deviations, if any, occurring 
at small k. The table layout emphasizes the qualitative cor- 
respondence between case (k.SS) and case {l/k,PS), ex- 
plained in Ref. 25. 

Some trends in the data are explained by appeaUng to the 
variational expression of the elastic energy W: 



W{e;k;f)^ inf 



(30) 



where JC{e) = {e; = {diUj + djUi)/2, (e) = e} is the set 
of admissible strain fields. E.g., for an incompressible mate- 
rial under SS loading (12 8 at and ( [30l l imply: 



e £i 



.^)W + (lA)(ei)«. 



(31) 



Hence the standard deviation S'_D(^'(e||) of the parallel com- 
ponent of the strain is essentially finite, since (e||)''^^ is, in 
agreement with the analytical expressions of the SDs in the 
next section to which we refer the reader for this discussion. 
Consider now another strain field e', solution for an anisotropy 
ratio k' > k. Using it as a trial field for problem dSOb with k 
provides one inequality. Duplicating the argument with k and 
k', and e, e' interchanged, yields after some easy algebra in- 
volving dSTT i: 

^^'""^''^ (Jao-Iw 

which entails ( |29] l for k' k. Thus, the standard deviation 
S'D^^-' {e±) of the transverse (PS) component of the strain field 
increases with fc at / fixed, consistently with Table Hill More- 
over, using (|32] | and the equality 9fc(£D'^-' = — fc9fc(£|)'"'^' 
[from ( |29] l1 shows that under SS loading (ey)^^^ is a decreas- 
ing function of fc. These considerations hold for any fixed 
microstructure. 

Analyzing EFT calculations at / = 0.1 for various values 
of fc in log-log plots (not shown), we observe that (for this /) 
the SDs behave as powers of k with numerical exponents close 
to 1/4 or 3/4: e.g., under SS loading, SD^^'>{e±) decays as 
/c3/4 when fc ^ 0, and blows up as fc^/^ when k oo; mean- 
while, S'L>(i)((T^) = SD^^^ei_)/k - fc-i/4 as fc ^ and 
^ /j;-3/4 as fc — > OO. The "soft" case fc is in agreement 
with the dilute analytical expressions ( I33bb and ( I33cb below, 
which indicates that the computed systems remained in the di- 
lute regime / ^ fr{k) ^ fc. On the other hand, the "hard" 
case fc ^ oo where strong strain localization takes place (see 
map C in Table |I| is consistent with ( |33b| ) and ( |33c| i only if 
we replace / by fr{k) ~ fc~^ in these expressions. Thus, 



here, SD^^^ {e±_) blows up [see ( I34dl i1. but behaves as though 
the system remained in the cross-over regime. This informa- 
tion, extracted numerically, is not contained in the expressions 
(I35dl ). ( I38bl ). for which we could only produce limiting values. 

Actually, in the limit fc oo, infinite SDs in the transverse 
component of the strain result from its concentration as Dirac 
lines (see Sec. 1111 Bl i. and are linked to discontinuities (jumps) 
in its parallel component. This results in a deformation pat- 
tern by a "rigid block sliding" mechanism, the "rigid blocks" 
being here connected parts of matter separated by discontinu- 
ity lines. This block-sliding effect only takes place provided 
that the strain jump lines have "percolated". Below "per- 
colation", sliding is impossible in a linear material and the 
transverse strain fluctuations described by (I33bl ) strongly in- 
crease with fc as SD^^^ (ej^) k^/"^. On the contrary, beyond 
"percolation", sliding takes place and SD'^^^e^^ ~ fc^/'' in- 
creases in a weaker way, since sliding makes deformation eas- 
ier Analogous properties are found under PS loading, pro- 
vided that fc is replaced by 1/fc: e.g., SD'^'^^ {ei_) is a decreas- 
ing function of fc and blows up in the hard loading mode as 
~ fc^^/"* when fc — !■ 0. 

Table |IV] shows numerical results for the strain and stress 
field averages and SDs, plotted for various anisotropy ratios 
fc = 0.01, 1 and 100. EMA estimates are provided for com- 
parisons. Except when SDs blow up at strong anisotropy, the 
EMA estimates are in good agreements with EFT results, for 
porosities up to / '--^ 0.4. It is worth observing that, in situ- 
ations of high anisotropy ratios, a change in the structure of 
the strip patterns in the material coincides with a change of 
concavity of the standard deviations SD^^'' ) of the parallel 
component of the strain field. For instance, when SS loading 
is applied at fc ^ 1, the bands cover the whole medium at 
/ w n/8 w 0.4. Around this value, the quantity 5'_D(^)(e||) 
changes from a concave to a convex function of /. Such a 
change also occurs at / w 7r/8 for the same SDs when fc ^ 1 
and PS loading is applied. EFT field maps then indicate that 
the structure of the strain pattern also undergoes an abrupt 
change at this point (with the appearance of thinner strips link- 
ing closest neighboring voids - not shown). 



2. Standard deviations in the dilute limit 

For completeness, we include the standard deviations in the 
dilute limit (of relevance to nonlinear EMAs) computed from 
the HS-EMA. 

For SS loading, we find for / <C fr{k) (low anisotropy or 
low porosity), 

SD{e{) - ^^(ctii) - /l/2fcl/^ (33a) 
SD{e^) ^ f^'^e'\ (33b) 
SD{a^) - /i/2fc-i/4^ (33c) 

where the /^/^ proportionality of the SDs goes along "classi- 
cal" 0(/) dilute corrections to the effective moduli. In limits 
of infinite anisotropy, where necessarily / 3> fr{k) = 0, the 
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TABLE III: 55 and PS loadings. Comparisons between EMA estimates (solid lines) and EFT results (dots) at porosity / = 0.1, for averages 
of the strain along the loading direction in each phase, and standard deviations (SD) of stress and strain components in the matrix, vs. matrix 
anisotropy ratio k — X/jj,. Strains and stresses are normalized by the appropriate macroscopic component in the loading direction (macroscopic 
strain e = (e), or stress a — (cr)). SDs in the voids are irrelevant. 
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HS estimates provide: 



SD{e\i) = SD{(7\i 



4^2 
7r3/4V3 



^3/4 



f 64 V2 V3^3/4- 



y37r9/4V3 4\/2 
S'L'(£±) = 0, SD{ai_) = oo, fc ^ 0, 



while exact results readi^ 



k~.0, (34a) 
(34b) 
fc^oo,(34c) 



29/4 

SD{e±) = oo, S'i:'(f7_L) = 0, k^oo, (34d) 



1 fc 0, (35a) 



SD{ei_) = 0, S'i:'(CT±) = oo, /c ^ 0, (35b) 

S'i:'(£||) = 5L>(a||) = (2//7r)^/'' , fc ^ oo, (35c) 

SD{ei_) = oo, S'i:»(CTx) =0, /c ^ oo. (35d) 

For PS loading, the HS-EMA gives for / < fr{k): 

SD{ei\) - 5L>((7||) - /^/^fc^^/", (36a) 



SD{ai_) ^ f^'^k^'\ 



(36b) 
(36c) 



TABLE IV: Average shear strains £e in the matrix (a = 1) and in voids (a = 2) vs. porosity /. Comparisons between EMA estimates (solid lines), EFT results (dots) and exact analytical 
results at = 0, oo(dash-dotted lines) for particular values of the matrix anisotropy ratio k = A//tt, in pure shear (PS) and simple shear (SS) loadings. The normalization is the same as for 
Table 111 (see legend). 
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For infinite anisotropy where / ^ fr (k) = 0, the HS esti- 
mates provide: 

SD{e^^) = SD{<m) = ^^^/i/4^ fc ^ 0, (37a) 
S'i:'(e^) = cx), S'L»(cr^) 0, k^O, (37b) 

29/4 



S'i:>(e||) = SDia«) 



+ 



^/37r3/4 



3/4 



fc^oo, (37c) 



33/2,^9/4 29/4 

SD{si_)^0, SD{(7^)^oo, k-^oo, (37d) 

while exact results arei^ 

S'D(eii) = S'i:'(CT||) = (//7r)^/^ /c -> 0, (38a) 
S'D(e_L)=oo, S'D(cr^) = 0, /c 0, (38b) 

99/4 

+ ,38c) 



21/48 \t: ' t:^ ^ 
SD{ej_) = 0, SD{aj_) ^ oo, /c -> oo. (38d) 

Thus, the HS-EMA correctly reproduces the scaling behavior 
of the SDs in all cases. 



VI. CONCLUDING DISCUSSION 

To summarize, we compared the results of the HS-EMA to 
EFT calculations, and showed that the agreement is excellent, 
even in the non-trivial case of localizing behavior, as far as 
effective moduli and averaged fields are concerned, and pro- 
vided that the void concentration lies below 0.3. This result 
is relevant to the study of non-linear effective-medium tech- 
niques: the latter involving both an anisotropic EMA, and 
a specific self-consistent linearization procedure (which de- 
termines the effective anisotropy of the former), the present 
study shows that should strong deviations between EET and 
EMA results in nonlinear media be observed (in the similar 
set-up of a periodic voided medium, and in similar condi- 
tions of porosity and of effective anisotropy), they ought be 
attributed to the linearization procedure rather than to the un- 
derlying liner EMA, even in limits of high effective anisotropy 
(determined by the field fluctuations in the nonlinear theory). 
Also, the present work provides a useful independent confir- 
mation of the involved analytical analysis of Ref. 25. 

As a by-product of the study, of possible practical applica- 
tions, we showed that by combining a regular lattice of voids 
(which makes the structure lighter) and an anisotropic matrix, 
properties could be tuned so as to make the overall medium 
elastically isotropic in plane strain. 

We also studied analytically the lattice sums which underly 
the EMA approach, and showed that they possess a scaling 



property which, in the dilute limit of small porosity and at high 
(but finite) anisotropy, allows for a cross-over between regular 
and singular porosity dependence of the effective medium. A 
length scale ^ was associated to this scaling, and interpreted 
as an effective heterogeneity size. It mathematically diverges 
in the limit of infinite anisotropy. However, its physically as- 
sociated counterpart being constrained by the finite size of the 
cell in the periodic medium, cross-over occurs when the ef- 
fective heterogeneities "percolate", i.e. when ^ is trivially of 
order one. This corresponds to a strongly correlated regime of 
strain localization bands spanning the system. 

Actually, scaling properties of lattice sums similar to 
the one considered here, have already been pointed out by 
Barber,*^ elaborating on Hall's work,''^ in a purely mathemat- 
ical context (in particular, no explanation in terms of length 
scales was given). Here, we make a connection between this 
phenomenon and strain localization in anisotropic elastic me- 
dia. Barber's paper also provides a means to compute the 
cross-over function. However, our lattice sums lead to tech- 
nical difficulties which preclude the straightforward obtention 
of a similar result, and we leave this issue for future work. 

Moreover, we found that under loading in a "hard" direction 
of the anisotropic medium, standard deviations of the trans- 
verse strain component blow up continuously as a power of 
the anisotropy ratio, as though the system remained in a cross- 
over regime. This absence of finite threshold for diverging 
fluctuations, and the above-described behavior, suggest the 
existence of a special type of continuous phase transition, of 
infinite order (called a "weak phase transition" by Hall,^^), 
here obviously without symmetry breaking. The presence 
of logarithmic terms in k (identified numerically in the Ap- 
pendix) also hints in this direction, since logarithmic correc- 
tions to scaling constitute a hallmark of transitions of infinite 
order. However, a random version of the system should be 
investigated before reaching definite conclusions. 

Finally, it was observed in Ref. 25 that the singular effec- 
tive moduli in the limit of infinite anisotropy are directly con- 
nected to the hyperbolic character of the governing equations 
in this limit. The very existence of a cross-over shows with- 
out ambiguity that the problem, although elliptic from a strict 
mathematical point of view, presents a "quasi-hyperbolic" 
character at short distances for high but finite anisotropy. 
This observation may be of relevance to theoretical investiga- 
tions of granular materials, for which a model with a similar 
anisotropic matrix has been proposed^ 
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APPENDIX: ASYMPTOTICS OF LATTICE SUMS AND 
DILUTE LIMITS 

We extract the dilute expansions of the lattice sums S'a,^ 
when fc — > 0, oo as follows. With the notation 



A{a,x) = [2Jl(27^aa;)/x]^ 
write the lattice sums, with p = {p^ + Py)^^"^, as 



(A.l) 



Sx ^ I J2 ^("^P) {Pl-Plf/D{pl (A.2) 
^ -J2 Aia,p)Aplpl/D{p), 



(A.3) 



Py>l 



where -D(p) = 4p^p^ + /c(p^ — Py) ■ """^^ above expres- 
sions explicitly acknowledge the fact that the principal diago- 
nal px = py does not contribute to S\, and that the Cartesian 
axis Px — does no contribute to 5^. Introducing the lattice 
sums 

5i(a) = - VA(a,n), ^2(0) = - E ^(«'^')' ^^-4) 

n>l Px>0 

p„>l 

provides for fc — > cx3, 0, 

TT ^ — ' / 

Pm>0 
Py >^iPy "^Px 

= - E ^(«'^') = ^2(0) - 5i(a). (A.5) 



In 5^^°° (resp. 5*^"*") the contribution of the principal diag- 
onal (resp. Cartesian axis px = 0) has been subtracted from 
£'2(0). The factors \/2 and 1/2 result from n being replaced 
by \/2n on the main diagonal. 

One privileged tool for exact asymptotic expansions is the 
Mellin transformi^ The transform and its inverse are defined 
by: 



M[/(x);z] = 




1 

2iTT 



dzx-'M[f{x);z], 



where c lies within the analyticity strip (parallel to the imagi- 
nary axis) of M[f{x); z] in the z-plane. Shifting the inversion 
contour to the left encircles the poles on the negative z axis 
and provides the asymptotic series expansion around x — 
in positive powers of x. Conversely, shifting the contour to 
the right provides the asymptotic expansion near a; = cx) in 
negative powers of x. The Mellin transform provides, for 
< c < 3, 



A{a,p) = ^ 



^+*°° dz 1 2 r(z/2)r(3/2-z/2) 
_,^(27ra)-^7^r(2-z/2)r(3-z/2) 



Next appealing to the definition of the Zeta function for z > 1, 

C{z) = X]„>i 1/"^' ™d to Hai-dy's lattice sui"^*'-^^ 



E 



'0 {pI+pIY^' 



= az/2)f3iz/2), {z > 2) (A.6) 



where P{z) = E„>o(-l)"/(2"- + 1) is the Dirichlet (or 
Catalan) function^^ and interchanging the lattice sums and 
the contour integral yields: 



^i(a) 



-I dz r(z/2)r(3/2-z/2) 

:4 [^) 



^'/Vci-zoo (2™)- " ^ ^ r(2 - z/2)r(3 - z/2) ' 

(A.7a) 



r-5/2 



:Ciz/2)(3{z/2) 



c2-ioo(27ra)" 

r(z/2)r(3/2- z/2) 
^r(2-z/2)r(3-z/2)' 



(A.7b) 



where 1 < ci < 3 and where 2 < C2 < 3 as the result of 
the above restrictions. The following properties hold: r(z) 
has simple poles at negative integers z = — fc < and has no 
zeros; ({z) has only one simple pole at z = 1, and has (so- 
called "trivial") zeros at even, nonzero, negative integers; /3(z) 
has no poles and has zeros at odd negative integers. Then, by 
shifting the contour to the left in both integrals, only the poles 
z = 1, 0, contribute to 5*1, and only the poles at z — 2 and 
z = contribute to 82- Eventually we obtain: 



^i(a) 



32 
— ( 

'Stt 



2'Ka^ 



5*2 (a) = 1 — na^ 



(A.8a) 
(A.8b) 



The polynomial form of these expressions indicates that they 
are actually exact, since for such functions the asymptotic ex- 
pansion coincides with the function itself. 



kSJhf) 




FIG. 7: Quantity kS\{k,a) vs. void concentration / — ira'^ as 
k —f 00: comparison between l |A.9ct , dashed line, and a numerical 
calculation of the lattice sum, solid line, for k = 10®. Breakdown of 
lAJcll occurs at / = 7r/8 ~ 0.39. 
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The validity of ( lA.Sal i is linked to the decay of 

r(z/2)r(3/2- z/2) 



r(2-z/2)r(3-z/2) 



as Sftez 



We have in the limit /(z)- s-^\C{z)r{z/2)C{z) /T{-z/2)\. 
Owing to the reflection formula 

^-^/2r(z/2)c(z) - ^-(i-^)/2r((i - z)/2)c(i - z), 

|r(z/2)C(z)/r(-z/2)| - 7r^|C(-z)| n\ Hence the con- 
tribution of the integration line in the limit ci —oo is neg- 
ligible only if s < TT. In terms of a, this amounts to a < 1/2. 
The breakdown of the obtained expressions thus corresponds 
to the close-packing limit a ~ 1/2. A similar reasoning us- 
ing the corresponding reflection formula for (3{z) provides the 
same range of validity for 82- 

Combining jA.Sh ( IA.8I ) and ( fT9l ) then results in 



kS 



fe^O 



cfe^O 

-5^ 



32 „ 2 
— a — ZTTa , 

. 32 

1 a + TTC 

Stt 



Qk — >oo 



1 



16^2 



I6V2 



(a < 1/2) 
\ (a<l/2) 

)■ 



27ra^, (a < 



2\/2 



(A.9a) 
(A.9b) 

(A.9c) 

(A.9d) 



The restrictions a < 1/2 and a < l/(2\/2) correspond to 
critical concentrations / = 7r/4and/ — tt/S. At these points, 
either the voids percolate (k = 0, PS or SS and k — cx), SS) or 
the shear bands undertake a configurational change (k = 00, 
PS)i^ An illustration of the breakdown of expression ( |A.9cl ) 
is provided by Fig.|7l 



We could not compute analytically the leading corrections 
in k to these sums. However, by carefully analyzing brute 
force numerical computations of the sums for k down to 
10^^ or up to 10^ for / = 0.1, we found that corrections 
to ( |A.9ab , ( |A.9b| ), ( |A.9ct . ( |A.9d| ) are of the form 0(A:logfc), 
0{—k log fc), 0(^k~^ log k), 0[~k~^ log k), respectively. 



As a final remark, we emphasize that only (divergent) 
asymptotic series for 5*1(0), 5*2(0) at a > 1/2 can be ob- 
tained: then, the integrand in both contour integrals blows up 
as [z / {4:TTae)Y where e is Euler's constant. These asymptotic 
expansions are easily extracted. We do not provide them here 
since the region a > 1/2 cannot be examined without appeal- 
ing to additional investigation procedures (e.g.. Fade approxi- 
mants) which lie outside the scope of this paper. 



